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ABSTRACT 


Atmospheric  pressure  and  wind  stress  fluctuations  are  strongly 
coupled  to  sea  level  fluctuations  along  the  Carolina  coasts  at 
periods  of  2.5  to  3.5  days.  Sea  level  fluctuations  in  this  band 
exhibit  high  coherence  over  a horizontal  separation  exceeding  500  km. 

Phase  difference  calculations  indicate  soutliward  propagation  of  the 
sea  level  fluctuations  from  Beaufort  to  Wilmington,  North  Carolina; 

1 

the  data  considered  are  insufficient  to  conclusively  determine 
propagation  direction  south  of  Wilmington.  The  2.5  - 3.5  day  period 
sea  level  fluctuations  are  consistent  with  a theoretically  expected 
first  mode,  barotropic  continental  shelf  wave.  It  is  concluded  that 
continental  shelf  waves  forced  by  the  atmospliere  contribute  to  the 
shelf  and  slope  water  circulation  off  North  Carolina. 


i 
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1.  Introduction 


It  is  commonly  observed  that  the  movements  of  coastal  and  continental 
shelf  waters  are  related  to  atmospheric  variables  such  as  the  local  wind. 

Over  the  last  decade  or  so  it  has  gradually  been  realized  that  a class  ; 

of  low  (subinertial ) frequency  vjave  motions  known  as  continental  shelf  i 

waves  (CSW's)  provide  a physical  mechanism  bv  means  of  whicli  the  energy  of 

atmospheric  fluctuations  can  be  focused  into  energetic  continental  shelf 

and  slope  water  motions.  Lov;  frequency  sea  level  fluctuations  associated 

with  these  motions  are  often  observed  to  be  coherent  over  hundreds  of 

kilometers  in  the  alongshore  direction.  Examples  of  the  existence  of 

oceanic  fluctuations  in  this  class,  detected  as  low  frequency  sea  level 

fluctuations,  are  given  for  the  Australian  coasts  by  Hamon  (1962),  for 

the  western  U.S.  coast  by  Mooers  and  Smith  (1968),  for  the  fiortli  Carolina 

coast  by  Mysak  and  Hamon  (1969) , and  for  the  east  Florida  coast  by 

Brooks  and  Mooers  (1977).  The  details  of  atmosphere-ocean  coupling 

processes  and  tiie  degree  of  their  Importance  to  the  total  circulation 

over  continental  shelves  are  still  a m.atter  of  active  debate,  but  there 

is  little  doubt  that  CSW's  universally  account  for  at  least  part  of  tlie 

low  frequency  variability  of  coastal  and  shelf  waters.  The  currently 

mounting  pressure  for  prudent  coastal  and  shelf  water  management  strategies 

calls  for  an  improved  understanding  of  how  CSW's  are  forced  by  the 

atmosphere  and  how  they  influence  the  circulation  patterns  over  the  shelf. 

This  study  examines  linkages  bctv'een  atmospheric  forces  and  low  frequency 
sea  level  fluctuations  along  the  Carolina  coasts,  taking  account  of  the 
coherent  relationship  between  surface  atmosniieric  pressure  and  winds. 


The  results  of  the  study  are  discussed  in  light  of  CSW  parameters  sucli 
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as  period  and  wavelength  expected  for  the  North  Carolina  shelf  region 
(Brooks,  1976a;  hereafter  referred  to  as  R76) . 

The  Carolina  shelf  is  the  northeastward  arm  of  the  South  Atlantic 
Bight,  Fig.  1.  The  Cape  Hatteras  region  of  the  North  Carolina  shelf 
may  be  a singularly  important  zone  of  CSW  generation,  since  the  departure 
of  the  Gulf  Stream  from  the  coast  there  represents  an  abrupt  change  in 
the  continental  shelf  waveguide  characteristics. 

In  the  presence  of  the  background  relative  vortlcity  of  the  Gulf 
Stream  off  Cape  Fear,  low  mode  stable  barotropic  CSW's  are  expected  to 
propagate  southward,  in  opposition  to  the  advectiv'e  tendency  of  the 
Stream,  (B76)  . The  zero  group  speed,  first  mode  CS\>.'  (likely  to  be 
forced  by  rapidly  traveling  atmopshcric  cold  fronts)  has  a period  of 
about  three  days  and  an  alongshore  wavelength  of  about  20  km;  its 
southward  phase  speed  is  thus  140  km  day  ^ (160  cm  s ^).  The  period  and 
wavelength  were  found  to  be  quite  sensitive  to  the  details  of  the  current 
profile,  because  the  upstream,  propagation  speed  of  the  wave  is  of  the 
same  order  as  tiie  downstream  speed  of  the  current.  The  predicted  along- 
shore wave  velocity  component  (v)  has  a maximum  value  near  the  coast, 
a node  near  the  sliclf  break  (approximately  100  km  offshore  of  Cape  Fear), 
and  decays  rapidly  seaward  of  the  shelf  break.  The  predicted  cross 
shelf  wave  velocity  component  (u)  has  a maximum  value  near  the  shelf 
break,  and  decreases  seaward  and  slioreward  of  the  shelf  break. 

The  right-handed,  "natural"  coordinate  system  used  in  B76  and  in 
this  study  has  been  rotated  clockwise  55.6  degrees,  such  that  the  x 
axis  is  appro.ximatel  V orthogonal  to  the  isobaths  off  Cape  Fear,  Fig.  2, 
line  A-B;  this  is  the  line  along  which  Richardson,  Schmitz,  and  Niiler 


5 


f 

(1969)  performed  a two-month  dropsonde  measurement  of  the  Gulf  Stream's 
velocity  structure.  The  coordinate  origin  is  point  A. 
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2 . Data  Processlns 

The  data  analyzed  are  from  coastal  sea  level  (tide  gauge)  and 
meteorological  stations  located  in  Figs.  1 and  2.  Hourly  values  of  sea 
level  height  for  1974  were  obtained  from  National  Ocean  Survey,  NOAA , 
Rockville,  Maryland,  for  stations  at  Beaufort  (BFT)  and  Wilmington  (WIL) , 
North  Carolina;  and  at  Charleston  (CHS),  South  Carolina.  Three-hourly 
values  of  surface  wind  speed,  wind  direction,  and  atmospheric  pressure 
for  1974  were  obtained  from  the  National  Climatic  Center,  NOAA,  Asheville, 
North  Carolina  for  stations  at  Cape  Hattcras  (HAT)^,  North  Carolina,  WIT, 
and  CHS. 

Tlie  raw  data  set  chosen  for  this  study  is  small  subset  of  the  total 

amount  of  coastal  atmospheric  and  tide  gauge  data  available  in  the  South 

Atlantic  Bight  from  the  year  1974.  The  stations  BFT  and  WIL  were  selected 

to  compare  with  a similar  study  of  earlier  data  from  stations  near  those 

locations,  Hysak  and  llamon  (1969).  CHS  was  Included  as  a third  station 

2 

because  of  its  availability  at  the  beginning  of  the  study  , and  to  increase 
the  maximum  available  alongshore  separation  dlst.ance  (BFT-CHS)  to  about 
500  km. 

3 

The  sea  level  data  were  low  pass  filtered  to  attenuate  the  daily 

and  semidaily  tides  and  Inertial  fluctuations  (the  inertial  period  at 

Cape  Fear  is  21.5  hours).  The  envelope  of  the  "iwo-day  low  pass"  filter 
energy  response  function  is  shown  in  Fig.  3.  .Attenuation  at  diurnal  and 
iiigher  frequencies  is  greater  than  10^'.  After  filtering,  the  sea  level 
data  were  resampled  at  8-hour  intervals. 

^Atmospheric  data  were  not  available  at  Beaufort. 

2jt  was  subsequently  learned  that  a tide  gauge  at  Avon,  N.  C.  (near 
Cape  Hatteras)  was  operational  during  1974;  the  data  from  Avon  will  be 
included  in  a subsequent  study. 

3a  I.anczos  filter  taper  v;ith  241  weights  (10  days  in  length)  was  used.  Thus 

five  days  (120  data  points)  are  lost  from  the  beginning  and  from  tlie  end  of 

the  raw  time  series.  The  low  nass  filtered  time  series  start  at  0000  hrs 


Fig.  3 Low  pass  filter  energy  response  envelope  function.  Attenuation 
at  1 CPD  is  10^. 


Three-hourly  wind  stress  vector  components  in  the  rotated  coordinate 
system  were  computed  from  the  raw  wind  speed  and  direction  data,  with  the 
positive  vector  sense  in  the  direction  toward  which  the  wind  blows.  The 
stress  components  were  computed  using  a quadratic  drag  lav;  with  the  drag 

-3 

coefficient  = 1.5  x 10  . The  wind  stress  components  and  atmtispheric 

pressure  time  series  were  low  pass  filtered,  using  a filter  witli  the 
response  characteristics  shown  in  Fig.  3.  The  atmospheric  data  were  then 
subsanipled  at  9-hour  intervals  and  linearly  interpolated  to  8-hour  intervals 
to  be  commensurate  with  the  sea  level  data. 

The  analyses  were  conducted  for  two  subsets  of  the  filtered  1974  data. 

A 120-day  subset  (6  Jan  to  6 May)  was  chosen  to  reflect  late  winter-early 
I spring  atmospheric  forcing  conditions  (i.e.,  to  emphasize  the  relative 

importance  of  cold  front  passages  over  the  Carolina  shelf  during  that 
period,  compared  to  the  distinctively  different  character  of  summer 
t atmospheric  conditions).  A 292-day  subset,  6 Jan  to  25  Oct,  was  chosen 

^ to  permit  an  approximate  comparison  with  the  1953  analysis  period  used  by 

I 

Mysak  and  Hamon  (1969)  (2  Jan  to  4 Sept);  and  to  provide  coverage  of  a 
^ period  in  which  "winter-like"  and  "summer-like"  atmospheric  forcing 

conditions  were  given  roughly  equal  time  weights. 

f 

i Extensive  use  is  made  of  spectrum  analysis.  Except  where  noted 

( 

otherwise,  all  frequency  domain  calculations  were  made  by  Fourier  trans- 
forming tapered  correlation  functions.  The  effective  spectrum  bandwidtli 
is  0.0667  cycles  per  day  (CPD) , with  16  (39)  degrees  of  freedom  for  the 
120  (292)  day  subset.  "Coherence  squared"  has  been  abbreviated  to  "coherence" 
in  the  text;  all  confidence  intervals  and  significance  tests  are  referred 
to  the  95%  level. 
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3 . Time  Domain  Analysis 

The  filtered  atmospheric  pressure  fluctuations  show  a high  degree  of 

4 

visual  correlation  between  the  CHS,  WTL,  and  Hj\T  stations  (Fig.  4)  . The 
similarity  is  apparent  even  for  relatively  small  amplitude,  short  period 
fluctuations  (e.g.,  days  42-46),  indicating  a horizontal  coherence  length 
scale  of  the  pressure  field  much  greater  than  the  CUS-HAT  Separation 
distance  (about  500  km) . The  large  horizontal  coherence  distance  pre- 
sumably reflects  the  synoptic  meteorological  scales  (thousands  of  l:m) 
associated  with  mid-latitude  winter  atmospheric  disturbances. 

The  alongshore  wind  stress  component  (Fig.  5)  also  shows  considerable 

correlation  over  the  CHS-llAT  separation  distance,  but  variations  in 

intensity  and  structure  between  stations  are  more  apparent  than  in  the 

_2 

pressure.  The  largest  stress  component  (+2.2  dyne  cm  ) occurred  at  HAT 
at  day  54;  this  "event"  can  be  clearly  identified  at  WIL  and  CHS,  where 
it  was  progressively  less  intense.  A vector  representation,  or  "stick" 
diagram,  of  the  wind  stress  shows  the  clockwise  rotation  of  the  t^ind 
field  usually  associated  with  the  passage  of  a cold  front  over  the 
station  (Fig.  6).  Again,  there  is  much  similarity  between  the  stations, 
with  the  largest  intensity  occurring  at  HAT.  The  strongest  winds  tend 
to  align  with  the  coastal  topography,  as  indicated  by  near-vertical  wind 
stress  vectors.  Wind  stress  "reversals"  (defined  by  a sign  change  of 
the  alongshore  component),  most  clearly  noted  at  HAT,  occurred  on  time 
scales  cf  several  days  to  several  v;ecks. 

Sea.  level  responded  to  the  major  wind  stress  events  in  the  manner 
expected  from  coastal  upwcl.ling  and  downwelling  (Fig.  7).  For  example, 
the  wind  stress  events  with  a strong  nf>r theastward  component  that  occurred 


The  origin  of  all  time  domain  figures  is  0000  hrs  01  January  1974. 
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Fig.  4 Filtered  atmosplierlc  pressure  at  tliree  coastal  stations 
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Filtered  wind  stress  vectors  In  the  rotated  coordinate  system 
for  three  coastal  stations. 
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Fig.  7 Wind  stress  vectors  at  Cape  Hatteras  and  their  relation  to 

unadjusted  (solid  lines)  and  barometrically  adjusted  (dashed 
lines)  sea  level.  Static  adjustment  Is  accomplished  by 
adding  1.01  cm  from  sea  level  for  each  local  atmospheric 

pressure  increase  of  1 mb. 
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during  days  53-55  and  60-88  caused  a total  sea  level  set-down  of  about 
55  cm  at  BFT,  consistent  with  an  offshore  Ekman  transport  of  surface  waters 
and  coastal  upwelling.  The  abrupt  reversal  at  day  70  to  persistent  wind 
stress  with  a southwestward  component  during  days  70  - 74  (downwelling- 
favorable)  produced  a sea  level  set-up  of  about  30  cm  at  BFT.  Similar 
responses  of  different  magnitudes  occurred  at  CHS  and  WIL. 

A prominent  sea  level  oscillation  with  a several  day  period  occurred 
at  all  stations.  In  most  cases,  the  several  day  oscillation  is  enhanced 
in  the  barometrically  adjusted  sea  level  curves  (dashed  line),  indicating 
a selective  atmospheric  pressure  ocean  coupling  process  at  this  period. 
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4 . Frequency  domain  analysis 

In  this  section,  the  coherence  and  phase  between  various  pairs 
of  atmospheric  and  sea  level  variables  are  discussed.  First,  atmos- 
plierlc  pressure  and  wind  stress  are  examined;  second,  sea  level  is 
examined;  and  third,  some  linkages  between  the  atmosphere  and  sea 
level  are  examined.  All  calculations  in  this  section  were  carried 
out  for  only  one  input  function,  i.e.  no  correction  has  been  made 
for  coherent  relationships  with  other  possible  input  functions. 

The  large  horizontal  coherence  scale  of  atmospheric  pressure 
fluctuations  deduced  from  the  time  domain  results  is  confirmed  by  the 
high  coherence  (-0.9)  between  pressure  at  HAT  and  CHS,  separated 
approximately  500  km  (Fig.  8B) . The  pressure  fluctuations  are  nearly 
in  phase  for  all  frequencies  in  the  range  0 to  0.5  cycles  per  day 
(CPD) , indicating  that  synoptic  scale  pressure  disturbances  arrive  at 
HAT  and  CHS  at  nearly  the  same  time.  The  largest  phase  lag^  occurred 
at  a period  of  about  3 days,  with  HAT  lagging  CHS  by  several  hours. 

Both  wind  stress  components  at  WiL  were  significantly  coherent  with 
the  pressure  for  periods  of  2.5  to  10  days,  with  the  cross  shelf  (along 
shelf)  component  exliibiting  a quadrature  (antiphase)  tendency  with 
respect  to  the  pressure. 

The  coherence  between  unadjusted  sea  level  records  for  tlie  120  day 
period  shows  rather  sharp  peaks  at  periods  of  2.5  - 3.5  and  7-30  days 
(Fig.  9).  The  highest  coherence  (0.96)  occurred  at  a period  of  3 days  in 
the  calculation  lietween  BFT  and  WIL,  confirming  the  time  domain  obser- 
vation of  prominent  several  day  oscillations  at  those  stations.  Surpris- 
ingly, the  coherence  at  periods  longer  than  10  days  is  higlier  between 


^ In  all  frequency-domain  figures,  positive  phase  angles  mean  that  the 
f i rs t-named  vir i ah  1 e leads  the  second-named  variable. 


Fig.  8 Coherence  and  phase  for  atmospheric  pressure  and  wind  stress 
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BFT  and  CHS  (the  "extrema"  stations)  than  between  WIL  (the  "central" 
station)  and  BFT  or  CHS;  this  may  reflect  the  fact  that  the  WIL  tide 
gauge  is  located  a considerable  distance  upstream  from  the  mouth  of 
the  Cape  Fear  river. 

Clearly  resolved  coherence  peaks  are  also  apparent  in  the  2.5  - 
3.5  and  7-30  day  bands  when  the  multi-season,  292-day  unadjusted 
records  are  considered  (Fig.  10).  The  phase  relationships,  which  are 
similar  for  both  record  lengths,  generally  indicate  fluctuations  at 
BFT  leading  those  at  WIL.  On  the  other  hand,  the  fluctuations  at  CHS 
lead  those  at  WIL  and  BFT  in  the  coherent  bands,  possibly  reflecting 
a static  sea  level  oscillation  traveling  northeastward  along  the  coast 
in  phase  with  synoptic-scale  atmospheric  pressure  systems. 

When  the  sea  level  is  statically  adjusted  for  atmospheric 
pressure  (for  each  mb  of  pressure  increase  at  the  nearest  atmospheric 
station,  1.01  cm  is  adding  from  the  sea  level),  the  coherence 
between  stations  increases  significantly  in  the  spectral  gap  between 
the  two  peaks  noted  earlier  (Fig.  11,  120  day  record;  Fig.  12,  292  day 
record).  However,  the  highest  coherence  values  consistently  occur  in 
the  2.5  - 3.5  day  band,  even  after  the  static  correction  is  made.  In 
that  band,  the  phase  indicates  fluctuations  at  BFT  leading  those  at 
WIL  by  0.3  to  0.4  day.  The  phase  results  in  the  2.5  to  3.5  day  band 
indicate  CHS  leading  WIL  and  BFT  by  about  0.3  and  0.1  days,  respec- 
tively; the  lead  is  somewhat  smaller  than  in  the  cases  of  unadjusted  sea 
level . 

The  cross  shelf  and  the  along  shelf  wind  stress  components  are  both 
consistently  coherent  with  the  120-day  record  adjusted  sea  level 


fluctuations  in  tlie  2.5  - 3.5  day  band  (Fig.  13,  cross-shelf  component; 
Fig.  14,  along  shelf  component).  Tlie  coherence  is  generally  low  for 


Coherence  and  phase  for  unadjusted  292-day  sea  level  rec 
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Fig.  14  Coherence  and  phase  for  along  shelf  wind  stress  components  and 
barometrically  adjusted  120-day  sea  level  records. 
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both  components  at  periods  longer  than  about  10  days,  although  the 
case  for  along  shelf  wind  stress  at  HAT  (Fig.  14C)  may  be  an  exception 
due  to  tlie  relative  exposure  of  that  station  to  the  Mid-Atlantic  Bight. 

The  phase  relationsliips  generally  indicate  falling  sea  level  witli 
increasing  cross  shelf  and  along  shelf  Vv'lnd  stress;  consistent 
respectively  with  coastal  sea  level  set-down  due  to  offshore  winds, 
and  v^'lth  offshore  Ekman  transport  and  resultant  set-down  due  to  upwel  ling- 
favorable  nortlieastward  winds. 

The  emergence  of  the  2.5  - 3.5  day  period  band  as  an  important 
portion  of  the  atmosphere-ocean  coupling  spectrum  indicates  a selective 
forcing  of  the  zero  group  speed,  first  mode  CSW  for  the  Cape  Fear  area, 
wliich  has  a period  of  3.0  days  and  a wavelength  of  420  km  (B76) . This 
CWS  is  predicted  to  propagate  southward,  in  opposition  to  the  Cult 
■Stream,  which  should  produce  lagging  phases  of  sea  level  f 1 uctu.at  ions 
associated  with  the  wave  as  one  proceeds  southward  along  the  coast. 

Hie  phase  evidence  presented  indicates  soutliward  propagation  from  BFT 
to  U'll,  (the  implied  wavelength  ,at  the  3-day  period  is  800-1000  km,  a 
t’ ictor  of  two  or  more  larger  tli.an  predicted),  but  the  phase  evidence 
for  soutliward  propagation  is  at  best  ambiguous  between  CHS  and  tlie 
other  stations.  Addition.il  se.i  level  stations  (me  or  two)  located 
between  CHS  and  WII,  would  probahlv  resolve  the  amhiguitv.  .\n  addition.il 
uncertainty  concerning  direction  of  phasi'  propagation  is  introduced  by 
the  slieltered,  up-river  location  of  the  Is'll.  tide  gauge.  I’.art  oi  the 
observed  BFT  and  WII,  pli.ise  lag,  may  be  due  to  the  time  required  for  the 
disturbincp  to  p rop.ig.i to  up  the  river  to  the  gauge  location.  Ho\^’ever, 
Mysak  ami  Hamt'ii  ( I ) also  found  southward  phase  propagation  in  Onslow 
Bay,  confirming  the  results  of  this  studv,  using,  a gauge  located  .it 
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Southport,  near  the  river  mouth.  (The  Southport  tide  gauge  was  not 
in  operation  during  the  period  of  this  study). 

In  summary,  the  frequency  domain  analyses  indicate  a highly 
coherent,  selective  coupling  process  between  the  atmosphere  and  sea 
level  in  the  2.5  --  3.5  day  band.  The  sea  level  fluctuations  were 
also  coherent  in  the  7-30  day  band,  but  coupling  with  atmospheric 
wind  stresses  was  weaker  than  in  the  2.5  - 3.5  day  band.  The  hori- 
zontal coherence  scales  of  atmospheric  and  sea  level  fluctuations 
w'ere  much  larger  than  the  greatest  station  spacing  in  this  study  (-500  km). 
The  phase  calculations  suggest  soutiiward  CSW  propagation  from  BFT  to 
U'lL,  consistent  with  Mysak  and  Hamon ' s results,  but  the  data  considered 
are  too  sparse  to  confirm  or  contraindicate  southward  propagation  along 
the  coast  south  of  K'lL.  A study  using  more  closely  tide  gauges  is 
necessary  to  clearly  resolve  the  phase  propagation  direction. 
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5 . Multiple  Regression  Models 

It  Is  clear  from  the  preceedlng  time  and  frequency  domain  analyses 
that  sea  level  fluctuations  are  coupled  with  atmospheric  forces  In  pre- 
ferred frequency  bands.  Since  atmospheric  pressure  and  winds  may  be 
correlated  with  each  other  to  a significant  degree  in  the  same  frequency 
bands.  It  is  desirable  to  correct  the  transfer  function  between  sea  level 
and  a particular  atmospheric  forcing  variable  for  the  coherent  contribution 
of  all  other  included  forcing  variables.  Thus,  for  example,  the  transfer 
function  between  atmospheric  pressure  and  sea  level  (the  "barometric 
function")  can  be  corrected  linearly  in  each  frequency  band  to  remove  the 
effects  of  wind  stresses  that  are  coherent  with  the  pressure.  Failure  to 
recognize  correlation  between  atmospheric  forcing  variables  can  cultivate 
misleading  physical  models;  e.g.,  that  sea  level  fluctuations  in  a parti- 
cular frequency  band  are  entirely  forced  by  pressure  fluctuations,  when 
in  fact  significant  forcin'’  may  be  duo  to  a wint!  stress  component  which  is 
highly  colierent  with  the  pressure  in  that  hand. 

The  commonly  applied  static  barometric  sea  level  "correction" 

( c f . , Mysak  and  Hamon,  1969)  is  based  on  the  premise  that  at  sufficiently 
low  frequencies  sea  level  responds  as  an  "inverse  barometer"  to  changes 
in  pressure  (a  pressure  increase  of  1 mb  results  in  a static  sea  level 
depression  of  1.01  cm).  The  static  correction  usually  does  account  for  a 
portion  of  local  sea  level  f li.ictuations , but  a significant  nonbarometric 
residual  correlation  often  remains,  indicating  the  importance  of  dynamic 
as  well  as  static  processes  (Brooks  and  Mooers , 1977).  The  static  correction 
seems  to  bo  a iielpful  first  step  in  exiioslng  dynamic  relationships  between 
atmospheric  pressure  and  sea  level  fluctuations.  A more  detailed  examination 
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should  at  least  recognize  the  likely  coherence  between  synoptic-scale 
atmospheric  pressure  and  wind  fluctuations,  particularly  since  CSW's  can 
apparently  be  forced  by  either  (Mysak,  1967;  Adams  and  Buchwald , 1969). 

Sea  level  fluctuations  are  considered  to  be  jointly  forced  by 
atmospheric  pressure  and  wind  stress  in  the  model  considered  in  this 
section.  The  linkages  between  the  variables  are  examined  by  a multiple- 
input,  linear,  frequency-domain  regression  scheme  (Groves  and  Hannan, 

1968).  Briefly,  the  simultaneous  regression  of  an  "output"  function  Y 
on  several  "input"  functions  Is  carried  out  in  frequency  bands  by 
operations  performed  on  the  augmented  cross  spectral  matrix  containing 
all  combinations  of  Y and  The  ordinary  coherence  and  phase  between 

the  output  and  each  input  are  linearly  corrected  for  the  coherent  effects 
of  all  tlie  other  inputs  in  each  frequency  band  to  yield  the  multiple 
coherence.  The  transfer  function  between  the  output  and  each  input  is 
likewise  corrected  for  all  other  coherent  inputs.  The  calculations  were 
carried  out  with  16  (47)  degrees  of  freedom  for  the  120  (292)  day  period. 

The  computer  programs  used  to  perform  the  calculations  are  part  of  the 
:<CSL'  FESTSA  time  series  analysis  system  (Brooks,  1976b). 

Multiple  coherence  between  sea  level  at  BFT  and  WIL  and  atmospheric 
variables  is  shown  in  Fig.  15  for  120  and  292  dav  cases.  Each  graph  shows 
the  coherence  increase  as  pressure  (p) , offshore  wind  stress  (t^)  and 
alongshore  wind  stress  (t^)  are  cumulatively  included  in  the  calculation. 

The  horizontal  dastied  line  is  tiie  95’1  null  hypothesis  level.  The  vertical 
arrows  identify  the  frequencies  of  the  first  and  second  mode,  zero  group 
speed  CSW's  determined  for  the  Cape  Fear  section  (B76) . Figure  15  (A  and  B) 
also  shows  the  spectrum  density  for  Gape  Hatteras  atmospheric  pressure  for 


the  120  and  292  day  period,  respectively. 


Fig.  15  Multiple  coherence  between  sea  level  and  a)  atmospheric  pressure, 
b)  pressure  and  cross  shelf  wind  stress,  and  c)  pressure,  cross, 
and  along  shelf  wind  stress  for:  A)  BFT,  120-day  record; 

B)  BFT,  292-day  record;  C)  Wll,  120-day  record;  and  D)  WIL, 
292-day  record.  Arrows  mark  zero  group  speed  frequencies  for 
first  and  second  CSW  modes  (from  Brooks,  1976a).  Tlie  Cape 
Hatteras  pressure  spectrum  for  120-day  (A)  and  292-day  (B)  cases 
is  shown  by  the  heavy  dashed  line. 
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By  including  the  wind  stress  components,  the  coherence  is  increased 
substantially  at  both  stations  and  for  both  record  lengths.  Coherence 
peaks  are  prominent  and  significant  at  both  stations  in  period  bands  of 

2.5  to  3.5  days  and  7 to  30  days,  with  the  highest  coherence  value  (0.95) 
occurring  at  3.0  days  at  WIL  for  the  120  day  (winter)  record.  Tlie  2.5  to 

3.5  day  band  straddles  the  first  mode,  zero  group  speed  CSW  frequency, 
suggesting  that  CSW's  are  selectively  forced  in  that  band.  Barotropic 
CSW's  forced  in  the  7 to  30  day  band  should  exhibit  southward  phase  and 
group  speeds,  but  it  is  not  suggested  by  the  calculations  in  B76  that  there 
should  be  a selective  coupling  in  that  band.  The  7 to  30  day  peak  is 
stronger  relative  to  the  2.5  to  3.5  day  peak  when  only  the  winter  months 
are  included,  Fig.  15  A and  C,  probably  reflecting  a higher  degree  of 
organization  of  winter  atmospheric  forcing  patterns.  The  pressure  spectra 
also  reflect  this  difference,  with  the  winter  case  (Panel  A)  showing 
prominent  peaks  at  periods  of  about  6 days  and  3 days;  whereas  the  combined 
v,finter  and  summer  case  (Panel  B)  is  generally  "red,"  with  only  hints  of 
peaks  in  similar  bands.  The  6-day  winter  pressure  spectrum  peak  is 
probably  associated  with  atmospheric  cold  front  passages;  the  3 day  peak, 
which  is  less  energetic  by  a factor  of  three  (but  much  more  strongly  coupled 
to  sea  level),  is  probably  second  harmonic  evidence  of  the  impulsive  time 
scale  (~1  day)  of  the  pressure  and  wind  disturbances  usually  associated 
with  a cold  front  passage  over  a fixed  station  (Brooks  and  Mooers,  1977). 

The  WIl,  coherence  differs  substantially  from  BFT  only  in  the  period 
band  of  7 to  30  days.  In  this  band,  the  WIL  peak  is  less  well  defined  and 
indicates  relatively  less  accrual  of  coherence  as  the  wind  stress  components 
are  included.  The  smaller  WIL  coherence  at  periods  longer  than  7 davs  may 
occur  because  the  tide  gauge  Is  located  approximately  35  km  upstream  from 
the  mouth  of  the  Cape  Fear  river,  in  a slieltered  location  at  which  exposure 


to  fluctuations  occurring  on  the  shelf  offshore  nuiy  be  reduced. 

Superbaronet ric  (greater  in  magnitude  than  1.01  cm/mb)  sea  level 
responses  to  pressure  fluctuations  generally  occur  in  bands  of  high 
coherence  (Fig.  16).  The  response  is  generally  sub-barometric  in  bands 
of  low  coherence  when  the  effects  of  wind  stress  are  included.  The 
uncorrected  barometric  moduli  (X's)  at  K'lL  tend  to  overestimate  the 
corrected  moduli  for  periods  shorter  than  about  4 days;  a similar  result 
is  not  evident  at  BFT.  With  the  exception  of  the  winter  WIT  case,  the 
barometric  functions  generally  indicate  an  out-of-phase  relationship 
between  atmospheric  pressure  and  sea  level,  as  expected.  Phase  estimates 
are  shown  only  for  frequencies  at  which  the  95?  confidence  interval 
(which  is  proportional  to  the  product  of  the  coherence  and  barometric 
function  modulus)  is  less  than  360°. 

Mysak  and  Hamon  (1969)  show  uncorrected  barometric  functions,  b(f), 

for  Southport  (near  WIL)  and  Morehead  City  (near  BFT)  , for  the  years 

1953  and  1954.  Their  results,  presented  as  b(f)  trajectories  in  the 

complex  plane  with  frequency  as  parameter,  are  qualitatively  consistent 

with  the  uncorrected  moduli  shown  in  Fig.  16.  However,  when  the  coherent 

wind  stresses  are  included,  the  curious  rotary  phase  patterns  shown  by 

Mysak  and  Hamon  are  not  reproduced;  instead,  the  corrected  phases  quite 

consistently  indicate  the  expected  out-of-phase  relationship,  at  least 

in  bands  of  relatively  high  coherence.  The  moduli  most  strongly  modified 

2 

by  inclusion  of  the  wind  stresses  were  those  at  WIL  for  periods  shorter 
than  about  4 days.  The  largest  difference  between  corrected  and  uncorrected 


This  may  be  partially  due  to  .shallow  water  effects,  such  as  wind  set-up, 
in  the  Cape  Fear  River. 
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F’tg.  16  R.irotncl  r Ic  function  modulus  (corrected  for  coherent  wind  stress 
solid  line;  iincorrectcd  broken  line)  nnd  phase  for:  A)  BFT, 

120-day  record;  B)  BFT,  292-day  record;  C)  WTL,  120-day  record; 
and  0)  WIL,  292-day  record.  Arrows  mark  zero  group  speed 
frequencies  of  first  and  second  CSW  modes  (from  Brooks,  1976a). 
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moduli  occuiTed  at  a period  of  about  2 days,  consistent  with  Mysak  and 
fiamon's  b(f)  graphs  which  show  the  largest  (superbarometric)  moduli  for 
periods  shorter  than  2.5  days. 

The  corrected  superbarometric  peaks  in  Fig.  16  consistently  occur  in 
the  high  coherence  bands  of  Fig.  15,  further  suggesting  dynamic  sea  level 
responses  to  atmospheric  forcing  in  these  bands.  The  ordinary  coherence 
between  wind  stress  and  sea  level  (Figs.  13  and  14)  tends  to  show  peaks 
in  the  3 day  band,  but  not  at  periods  longer  than  10  days,  suggesting  that 
pressure  and  wind  stress  fluctuations  associated  with  cold  fronts  are 
jointly  most  effective  in  forcing  sea  level  fluctuations  in  the  2. 5-3. 5 
day  band . 

Sub-barometric  sea  level  responses  generally  occurred  at  periods  of 
4-6  days.  Those  are  difficult  to  explain,  particularly  in  view  of  the 
antiphase  nature  of  the  response,  for  one  might  expect  a barometric  back- 
ground function  with  superimposed  dynamic  responses  in  bands  '..'here  the 
sea  level  "admittance”  is  iiigh;  for  example,  at  C.SW  zero  group  speed 
frequencies.  it  is  curious  that  the  smallest  moduli  occur  near  the 
frequency  of  the  most  prominent  wnnter  pressure  spectrum  peak^  , which  also 
corresponds  closely  to  the  second  CSW  mode  zero  group  speed  frequency; 
tiiis  suggests  strongly  preferential  excitation  of  the  first  CSW  mode 
relative  to  the  second  mode. 

In  summary,  tb.e  multiple  regression  model  indicates  a rather  selective 
coupling  mechanism  between  the  atmospliere  and  sea  level  in  2. 5-3. 5 and 
7-30  day  bands.  The  2. 5-3. 5 day  band  corresponds  well  with  the  ‘irst 
barotropic  CSW  zero  group  speed  frequency:  the  winter  pressure  sju'ctriim 

^Mysal  and  H.amon's  smallest  b(f)  moduli  also  occurred  near  a period  of 
5 days,  indicating  some  measure  i)f  stationarity  of  the  process  over  two 
decaiies . 
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suggests  a 3-day  period  forced  by  the  second  harmonic  of  the  pressure 
field  CSW  response.  The  7-30  day  coupling  is  not  clearly  explained  as 
a stable,  barotropic  CSW  response.  Wlien  coherent  wind  stresses  are 
included,  the  barometric  function  generally  exhibits  the  expected  180 
degree  phase  relutionslilp  with  sea  level  fluctuations.  The  barometric 
moduli  also  suggest  the  occurrence  of  dynamic  sea  level  responses  in  the 
2. 5-3. 5 and  7-30  day  bands. 


< 
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6 . Summary  and  Conclusions 

The  most  highly  coherent  sea  level-atmospheric  coupling  occurred 
in  the  2.5  to  3.5  day  period  band,  consistent  with  an  atmospherically- 
forced,  southward  propagating  continental  shelf  wave  (CSW)  expected  to 
occur  in  this  band  (Brooks,  1976a).  The  phase  relationships  indicated 
soutliward  propagation  from  Beaufort  to  Wilmington  in  the  2.5  to  3.5  day 

^ band,  but  southward  propagation  from  Wilmington  to  Charleston  was  not 

I 

^ confirmed. 

! The  multiple  coherence  between  sea  level  and  atmospheric  variables 

i 

j also  indicates  strong  coupling  in  the  2.5  to  3.5  and  7 to  30  day  bands. 

The  wind  stress  components  made  the  least  contribution  to  multiple 
coherence  in  the  2.5  to  3.5  day  band  in  the  winter  season;  which  suggests 
that,  when  coherence  between  winds  and  atmospheric  pressure  are  accounted 


for,  the  pressure  fluctuations  provide  the  primary  forcing  of  sea  level 
fluctuations  in  this  band.  For  other  bands  in  the  winter  and  for  all 
hands  in  the  joint  winter-summer  season,  addition  of  the  wind  stress 
components  significantly  increased  the  multiple  coherence.  The  dominant 
vjinter  response  in  the  2.5  to  3.5  day  band  may  be  a CSW  forced  by  the 
second  harmonic  energv  of  impulsive  pressure  and  wind  stress  fluctuations 
associated  with  the  passage  of  cold  fronts  over  the  area.  The  sea  level 
response  to  the  fundamental  atmospheric  pressure  oscillation  (at  a period 
of  about  6 days)  was  considerably  smaller  than  at  the  3 day  period,  which 
indicates  a strongly  preferential  excitation  of  the  first  CSW  mode. 


compared  to  the  second  mode. 

The  atmospheric  pressure-to-sea  level  transfer  functions  (the 


barometric  function),  Tjhen  corrected  for  coherent  wind  stress  effects. 
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generally  indicated  the  occurrence  of  dynamic,  antiphase  sea  level 
responses  to  pressure  fluctuations  at  periods  of  2 to  4 days  and 
greater  than  10  days.  The  static  "barometer  factor"  (-1.01  cm/mb) 
only  crudely  approximates  the  barometric  function  when  averaged  over 
0 to  0.5  cycles  per  day,  and  then  only  if  the  coherent  wind  stresses 
are  included.  The  uncorrected  barometric  function  at  Wilmington 
substantially  over-estimated  the  corrected  one  for  periods  shorter  ttian 
about  4 days.  The  curious  "rotary"  barometric  function  phases  reported 
by  Mysak  and  Hamon  (1969)  are  not  reproduced  here  when  the  wind  stress 
components  are  included  in  the  calculations. 

In  conclusion,  these  results  indicate  that  atmospherically  forced 
CSW's  contribute  to  the  total  circulation  of  the  continental  shelf  and 
slope  waters  off  Morth  Carolina.  The  most  important  spatial  (temporal) 
scale  of  motion  to  be  expected  solely  as  a result  of  atmospherically 
forced  CSW's  is  about  400  km  (3  days);  other  scales  of  motion  arc  not 
excluded.  The  relative  importance  of  C . circulation  compared  to  that 
driven  by  or  associated  with  other  meihanisms  (such  as,  for  example, 
instabilities  of  the  Gulf  Stream)  remains  an  open  question  of  great 
interest . 
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List  of  Figure  Captions 

Fig.  1 Map  of  South  Atlantic  Bight,  showing  locations  of  tide  gauge 
and  atmospheric  station  data  used  in  the  study. 

Fig.  2 Onslow  Bay  continental  shelf  topography,  showing  coordinate 
system  orientation. 

Fig.  3 Low  pass  filter  energy  response  envelope  function.  Attenuation 
at  1 CPD  is  10^. 

Fig.  4 Filtered  atmospheric  pressure  at  three  coastal  stations. 

Fig.  5 The  filtered  along  shore  v;ind  stress  component  at  three  coastal 
stations . 

Fig.  6 Filtered  wind  stress  vectors  in  the  rotated  coordinate  system 
for  three  coastal  stations. 

Fig.  7 Wind  stress  vectors  at  Cape  Hatteras  and  their  relation  to 

unadjusted  (solid  lines)  and  barometrically  adjusted  (dashed 
lines)  sea  level.  Static  adjustment  is  accomplished  by 
subtracting  1.01  cm  from  sea  level  for  each  local  atmospheric 
pressure  increase  of  1 mb. 

Fig.  8 Coherence  and  phase  for  atmospheric  pressure  and  wind  stress 
components . 

Fig.  9 Coherence  and  phase  for  unadjusted  120-day  sea  level  records. 

Fig.  10  Coherence  and  phase  for  unadjusted  292-day  sea  level  records. 

Fig.  11  Coherence  and  phase  for  barometrically  adjusted  120-day  sea 

level  records. 

Fig.  12  Coherence  and  phase  for  barometrically  adjusted  292-day  sea 
level  records. 

Fig.  13  Coherence  and  phase  for  cross  shelf  wind  stress  components  and 
barometrically  adjusted  120-day  sea  Level  records. 

Fig.  14  Coherence  and  phase  for  along  shelf  wind  stress  components  and 
barometrically  adjusted  120-day  sea  level  records. 

Fig.  15  Multiple  coherence  between  sea  level  and  a)  atmospheric  pressure, 
b)  pressure  and  cross  shelf  wind  stress,  and  c)  pressure,  cross, 
and  along  shelf  wind  stress  for;  A)  BFT,  120-day  record; 

B)  BFT,  292-dav  record;  C)  Wll,  120-day  record;  and  D)  WIL, 
292-day  record.  Arrows  mark  zero  group  speed  frequencies  for 
first  and  second  C.BW  modes  (from  Brooks,  1 976a).  The  Cape 
Hatteras  pressure  spectrum  for  120-day  (A)  and  292-d.ay  (B)  cases 
is  shown  hv  the  heavy  dashed  line. 
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Fig.  16  Barometric  function  modulus  (corrected  for  coherent  wind  stress, 
solid  line;  uncorrected  broken  line)  and  phase  for:  A)  BFT , 

120-day  record;  B)  BFT,  292-day  record;  C)  WIL,  120-day  record; 
and  0)  WIL,  292-day  record.  Arrows  mark  zero  group  speed 
frequencies  of  first  and  second  CSW  modes  (from  Brooks,  1976a). 
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